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APPLICATIONS FOR PATENT 
Inventors: Joseph Nir, Baruch Shayevich and Hanoch Cohen 

5 

Title: COHERENT PROCESSING OF SATELLITE SIGNALS TO 

LOCATE A MOBILE UNIT 

10 FIELD AND BACKGROUND OF THE INVENTION 

The present invention relates to positioning systems for portable radio 
transceivers, and, more particularly, to a positioning system for a mobile unit of a 
wireless communications network, based on signals from a satellite network such as 
the Global Positioning System (GPS) and on signals from the base stations of the 
15 wireless communications network. 

The Global Positioning System is a system of low Earth orbit satellites that 
transmit standard signals that can be used to establish the location of a user equipped 
with a suitable GPS receiver. For civilian applications, the signals are "C/A codes" 
that are pseudorandom noise (PN) sequences transmitted at a chip rate of 1023 KHz 
20 and a repetition period of 1023 chips, so that each frame of the PN sequence includes 
1023 chips. Each satellite is assigned its own PN sequence, and the PN sequences of 
the various satellites are mutually orthogonal. 

Superposed on the PN sequences transmitted by a GPS satellite is a Satellite 
Data Message, transmitted every 30 seconds at a rate of 50 bits per second. Each bit 
25 (+1 or -1) of the Satellite Data Message is modulated onto the satellite's signal by 
multiplying 20 consecutive frames of the PN sequence by the bit. The first 900 bits of 
the Satellite Data Message includes the satellite ephemeris and time model for the 
respective satellite. The remaining 600 bits of the Satellite Data Message include a 
portion of the GPS almanac, which is a 15,000 bit block of coarse ephemeris and time 
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model data for the entire GPS system. In addition, bits 1-8, bits 301-308, bits 601- 
608, bits 901-908 and bits 1201-1208 of every Satellite Data Message are identical 8- 
bit (160 millisecond) headers that are invariant in time and that are identical for all the 
GPS satellites; and bits 31-60, bits 331-360, bits 631-660, bits 931-960 and bits 1231- 
5 1260 of every Satellite Data Message are 30-bit (600 millisecond) handover words 
that are time-variant (these handover words include representations of the time of the 
week), but that, like the headers, are identical for all the GPS satellites. 

Conventionally, a GPS receiver acquires and tracks the signals from at least 
four GPS satellites, by correlating the received signal with the satellites' respective 
10 PN sequences and locking on to the correlation peaks. Once the satellites are acquired 
and tracked, the GPS receiver decodes the ephemeris and time model, for each 
acquired satellite, from the respective Satellite Data Messages. These models include 
sufficient ephemeris data to enable the GPS receiver to compute the satellite's 
positions. The correlation peaks obtained during the continued tracking of the 
15 satellites provide measured times of arrival of these PN sequence frames. The 
differences between an arbitrary reference time and measured times of arrival, 
multiplied by the speed of light, are pseudoranges p from the satellites to the GPS 
receiver. Typically, the reference time is the time at which the satellites commenced 
transmission of their respective PN sequences, as measured by the GPS receiver clock, 
20 which in general is offset from the GPS system clock by an unknown time offset. A 
pseudorange p is related to the true range R of the respective satellite by p=R+c h , 
where the range offset c b is the time offset, 7b, of the GPS receiver relative to GPS 
system time, multiplied by the speed of light c: c b =T 0 c. From these pseudoranges, and 
from the known positions of the satellites as functions of time, the position of the GPS 
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receiver is calculated by triangulation. Pseudoranges to at least four satellites are 
needed to solve at least four simultaneous equations of the form 
\s-r\ = p-c b 

where s is the position vector of a satellite and r=(x,y,z) is the position vector of the 

5 GPS receiver, for the three unknown Cartesian coordinates x, y, z of the GPS receiver 
and for c b . The satellites are sufficiently far from the GPS receiver that these 
equations can be linearized in x, y and z with no loss of accuracy. 

Several methods are known for increasing the efficiency with which a GPS 
receiver establishes its position and for reducing the power requirements of a GPS 

10 receiver. Schuchman et al., in U. S. Patent No. 5,365,450, which is incorporated by 
reference for all purposes as if fully set forth herein, teach the provision of the 
Satellite Data Messages to a GPS receiver integrated into a mobile unit, such as a 
cellular telephone, of a wireless communications network such as a cellular telephone 
network, by transmitting the Satellite Data Messages to the mobile unit from a base 

15 station of the network via the control channel of the network. The Satellite Data 
Message is 30 seconds long, so even under ideal reception conditions with parallel 
processing of the signals from all the satellites in view, it necessarily takes more than 
30 seconds to get a GPS position fix. Prior knowledge of the Satellite Data Message 
reduces this time to under 10 seconds. 

20 Krasner, in U. S. Patent No. 5,663,734, which is incorporated by reference for 

all purposes as if fully set forth herein, teaches a GPS receiver for a mobile unit in 
which, as in the GPS receiver of Schuchman et al., the Satellite Data Message is 
obtained by a wireless link to a base station, but then, instead of processing GPS 
signals in real time, the GPS receiver stores up to one second's worth of signals (1000 
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PN sequence frames per satellite), along with the initial time of arrival of the signals, 
and processes the stored signals. Groups of 5 to 10 frames each are summed and 
correlated with the PN sequences of satellites expected to be in view, and the resulting 
correlation functions are added incoherently. The summation over up to 1000 frames 

5 boosts the signal to noise ratio by a corresponding amount, and the post facto 
processing requires much less power than real time processing. 

Duffet-Smith et al, in PCT Application WO 99/21028, which is incorporated 
by reference for all purposes as if fully set forth herein, teach a system and method for 
locating a mobile receiver ("remote unit") of a wireless communications network. 

10 The remote unit and a base unit at a fixed location both receive signals on a control 
channel from three or more base transceiver stations (BTSs), and both correlate 
invariant portions of the control channel signals. Compact descriptions of the 
correlation peaks, such as low order polynomial fits, are transmitted to a central 
processor by both the base unit and the remote unit. The central processor recovers, 

15 from the correlation peaks, the time offsets of the BTS transmissions to the base and 
remote units. These time offsets are the functional equivalents of GPS pseudoranges. 
Given the (fixed) locations of the BTSs and the base unit, the central processor 
computes the position of the remote unit, essentially by triangulation, and transmits 
that location to the remote unit. 

20 A GPS positioning system is predicated on the GPS receiver having a clear 

line of sight to at least four GPS satellites. This often is not the case in an urban 
environment. Urban environments commonly have cellular telephone networks 
installed, so in principle the system taught by Duffet-Smith et al. could be used for 
locating a mobile transceiver configured as a cellular telephone; but cellular 


telephones generally are in simultaneous contact with at most two base stations. 
Sheynblatt, in U. S. Patent No. 5,999,124, Camp, in PCT Application W099/61934, 
and Waiters et al., in U. S. Patent No. 5,982,324, all teach methods of using signals 
from both GPS satellites and terrestrial BTSs to determine the location of a mobile 

5 receiver. These prior art methods generally rely on triangulation, as described above, 
using signals received simultaneously from a combined total of four or more GPS 
satellites and BTSs. As noted above, usually, at most two BTSs are "in view" 
simultaneously; and in a sufficiently cluttered and noisy urban environment, useable 
signals may be available from only one or two GPS satellites at any one time. These 

10 prior art methods also rely on all the BTSs being mutually synchronized, at least with 
each other if not with the GPS system. Such synchronization is optional under the 
GSM cellular standard. There is thus a widely recognized need for, and it would be 
highly advantageous to have, a system and method for locating a mobile unit of a 
wireless communications network, for example a mobile cellular telephone in a 

15 cluttered and noisy urban environment, that derives the location of the mobile unit 
from signals received sequentially from low Earth orbit satellites such as GPS 
satellites and from base stations of the communications network, without requiring 
that the BTSs be mutually synchronized. 

20 SUMMARY OF THE INVENTION 

According to the present invention there is provided a method for determining 
a pseudorange, and a rate of change thereof, to a beacon that transmits a signal 
including a plurality of blocks, each block including a plurality of frames of a 
pseudonoise sequence, each block being multiplied by a bit of a data sequence, the 
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signal being shifted relative to a nominal frequency by a frequency shift, the method 
including the steps of: (a) receiving the signal; (b) digitizing the received signal, 
thereby producing a digitized signal including a plurality of bits; (c) arranging the 
digitized signal as columns of an input matrix that includes a plurality of rows, each 
5 column including consecutive the bits of the digitized signal that correspond to an 
integral number of frames of the pseudonoise sequence; and (d) performing a discrete 
orthogonal transform on each row of the input matrix, thereby producing a 
transformed matrix. 

According to the present invention there is provided a method for determining 
10 a pseudorange, and a rate of change thereof, to each of a plurality of beacons, each 
beacon transmitting a respective signal including a plurality of blocks, each block 
including a plurality of frames of a pseudonoise sequence, each block being multiplied 
by a bit of a data sequence, the pseudonoise sequences and the data sequences being 
uniquely associated with respective satellites, all the pseudonoise sequences being of 
15 equal length, the pseudonoise sequences being mutually orthogonal, the method 
including the steps of: (a) receiving the transmitted signals collectively as a received 
signal; (b) digitizing the received signal, thereby producing a digitized signal 
including a plurality of bits; and (c) for each beacon, multiplying groups of the bits of 
the digitized signal, that correspond to the blocks of the signal transmitted by the each 
20 beacon, by respective bits of the data sequence of the each beacon. 

According to the present invention there is provided a receiver for receiving a 
signal transmitted by a beacon, the signal including a plurality of frames of a 
pseudonoise sequence, including: (a) an antenna for receiving the transmitted signal; 
(b) a mechanism for digitizing the received signal to produce a digitized signal 
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including a plurality of bits; (c) a memory for storing the digitized signal as columns 
of a matrix that includes a plurality of rows, each column including consecutive the 
bits of the digitized signal that correspond to an integral number of frames of the 
pseudonoise sequence; and (d) a processor for performing a discrete orthogonal 
5 transform on each of the rows of the matrix. 

According to the present invention there is provided a locator system for 
locating a mobile unit, including: (a) at least one beacon having a respective 
pseudonoise sequence and a respective data sequence, each data sequence including a 
plurality of bits, each at least one beacon operative to transmit a respective transmitted 
10 signal, each respective transmitted signal including a plurality of blocks, each block 
including a plurality of frames of the respective pseudonoise sequence, each block 
being multiplied by a bit of the respective data sequence, (b) a reference unit 
including: (i) a reference unit receiver for: (A) receiving the at least one transmitted 
signal as a reference unit received signal, and (B) recovering the at least one data 
15 sequence from the respective reference unit received signal, and (ii) a transmitter for 
transmitting the at least one data sequence to the at least one mobile unit; and (c) in 
the mobile unit: (i) a first mobile unit receiver for receiving the at least one data 
sequence, and (ii) a second mobile unit receiver including: (A) an antenna for 
receiving the at least one transmitted signal collectively as a mobile unit received 
20 signal, (B) a mechanism for digitizing the mobile unit received signal to produce a 
digitized signal including a plurality of bits, (C) a memory for storing, for each of the 
at least one beacon, an instance of the digitized signal, and (D) a processor for 
multiplying each of the at least one instance of the digitized signal by the bits of the 
respective data sequence. 
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According to the present invention there is provided a method for determining 
a location of a receiver, including the steps of: (a) providing a plurality of beacons 
having respective pseudonoise sequences, all the pseudonoise sequences being of 
equal length; (b) transmitting, by each beacon, a respective signal including a plurality 

5 of frames of the respective pseudonoise sequence; (c) receiving the transmitted signals 
collectively as a received signal, by the receiver; (d) inferring, for each beacon, a 
pseudorange and a rate of change of the pseudorange; and (e) inferring the location of 
the receiver from the pseudoranges and from the rates of change of the pseudoranges. 

According to the present invention there is provided a method for determining 

10 a pseudorange, and a rate of change thereof, to a beacon that transmits a signal 
including a plurality of blocks, each block including a plurality of frames of a 
pseudonoise sequence, each block being multiplied by a bit of a data sequence, the 
signal being shifted relative to a nominal frequency by a Doppler frequency shift, the 
method including the steps of: (a) receiving the signal; (b) digitizing the received 

15 signal, thereby producing a digitized signal including a plurality of bits; and (c) 
applying a matched filter algorithm to the digitized signal to extract therefrom the 
pseudorange and the rate of change of the pseudorange, the matched filter algorithm 
including: (i) demodulating the digitized signal relative to the data sequence. 

According to the present invention there is provided a method for locating a 

20 mobile unit, including the steps of: (a) providing a plurality of asynchronous 
stationary beacons, each stationary beacon transmitting a respective stationary beacon 
signal; (b) providing a plurality of mobile beacons moving on known respective 
trajectories, each mobile beacon transmitting a respective mobile beacon signal, all the 
mobile beacons being mutually synchronized; (c) for each of at least one of the 
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stationary beacons: (i) receiving the respective stationary beacon signal, by the mobile 
unit, and (ii) inferring, from the respective stationary beacon signal, at least one 
parameter related to a range to the each stationary beacon; (d) for each of at least one 
of the mobile beacons: (i) receiving the respective mobile beacon signal, by the 

5 mobile unit, and (ii) inferring, from the respective mobile beacon signal, at least one 
parameter related to a range to the each mobile beacon; (e) synchronizing the mobile 
unit and the stationary beacons to the mobile beacons by steps including: (i) 
providing, for each stationary beacon, a reference unit for measuring time offsets 
between the mobile beacons and the each stationary beacon, (ii) measuring the time 

10 offsets, by the reference units, and (iii) transmitting at least one the time offset from a 
respective the reference unit to the mobile unit, the inferring of the parameters being 
based on the synchronization; and (f) inferring, from the parameters, a location of the 
mobile unit. 

According to the present invention there is provided a method for tracking a 
15 mobile unit, including the steps of: (a) providing a plurality of beacons, each beacon 
transmitting a respective ranging signal; (b) obtaining an initial estimate of a location 
of the mobile unit; (c) receiving, by the mobile unit, one of the ranging signals; (d) 
inferring, from the ranging signal, at least one parameter related to a range to the 
beacon that transmitted the ranging signal; and (e) obtaining an updated estimate of 
20 the location of the mobile unit, based on an initial state estimate that includes the 
initial estimate of the location and on an update measurement that includes the at least 
one parameter. 

According to the present invention there is provided a locator system for 
locating a mobile unit, including: (a) a plurality of asynchronous stationary beacons, 
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each stationary beacon being operative to transmit a respective stationary beacon 
signal; (b) a plurality of mobile beacons, each mobile beacon operative to transmit a 
respective mobile beacon signal while moving on a known respective trajectory, all 
the mobile beacons being mutually synchronized; and (c) in the mobile unit: (i) a 
5 mobile unit clock for producing a mobile unit timing signal; (ii) a first mobile unit 
receiver for receiving and digitizing the stationary beacon signals; (iii) a second 
mobile unit receiver for receiving and digitizing the mobile beacon signals, thereby 
producing corresponding digitized mobile beacon signal; (iv) a memory for storing the 
digitized mobile beacon signal; (v) a mobile unit processor for transforming each of 
10 the stationary beacon signals into a corresponding stationary beacon pseudorange, and 
for transforming each of the digitized mobile beacon signals into a corresponding 
mobile beacon pseudorange subsequent to the storing of the each digitized mobile 
beacon signal in the memory; and (d) at least one reference unit including: (i) a first 
reference receiver for receiving the mobile beacon signals and inferring therefrom a 
15 common start time of the mobile beacon signals, (ii) a second reference unit receiver 
for receiving at least one of the stationary beacon signals from at least one respective 
the stationary beacon, (iii) a reference unit clock for producing a reference unit timing 
signal, (iv) a reference unit processor for inferring, for each at least one respective 
stationary beacon, a timing offset between a start time of the at least one respective 
20 stationary beacon signal and the common start time of the mobile beacon signals; and 
(iii) a transmitter for transmitting the at least one timing offset; the first mobile unit 
receiver being operative to receive the at least one timing offset from each at least one 
reference unit transmitter; the pseudoranges being with reference to both the mobile 
unit timing signal and the at least one timing offset. 
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According to the present invention there is provided a locator system for 
locating a mobile unit, including: (a) a plurality of beacons, each beacon transmitting 
a respective ranging signal; and (b) in the mobile unit: (i) a mechanism for obtaining 
an initial estimate of a location of the mobile unit, (ii) a receiver for receiving the 
5 ranging signals, (iii) a processor for inferring, from one of the ranging signal, at least 
one parameter related to a range to the beacon that transmitted the ranging signal, and 
(iv) an updating mechanism for updating the location estimate, based on an update 
measurement that includes the at least one parameter. 

According to the present invention there is provided a method for tracking a 
10 mobile unit, including the steps of: (a) providing a plurality of mobile beacons, each 
mobile beacon transmitting a respective ranging signal; (b) providing a plurality of 
stationary beacons, each stationary beacon transmitting a respective ranging signal; (c) 
obtaining an initial estimate of a location of the mobile unit; (d) receiving, by the 
mobile unit, a first ranging signal from a respective beacon of one of the pluralities of 
15 beacons; (e) inferring, from the first ranging signal, at least one parameter related to a 
range to the beacon that transmitted the first ranging signal; (f) obtaining a first 
updated estimate of the location of the mobile unit, based on the initial location and 
on the at least one parameter that is related to the range to the beacon that transmitted 
the first ranging signal; (g) receiving, by the mobile unit, a second ranging signal from 
20 a respective beacon of another of the pluralities of beacons; (h) inferring, from the 
second ranging signal, at least one parameter related to a range to the beacon that 
transmitted the second ranging signal; and (i) obtaining a second updated estimate of 
the location of the mobile unit, based on the first updated location and on the at least 
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one parameter that is related to the range to the beacon that transmitted the second 
ranging signal. 

According to the present invention there is provided a locator system for 
locating a mobile unit, including: (a) a plurality of mobile beacons, each mobile 
5 beacon transmitting a respective ranging signal; (b) a plurality of stationary beacons, 
each stationary beacon transmitting a respective ranging signal; and (c) in the mobile 
unit: (i) a mechanism for obtaining an initial estimate of a location of the mobile unit, 
(ii) a receiver for receiving the ranging signals, (iii) a processor for inferring, from 
each of the ranging signals, at least one parameter related to a range to the beacon that 
10 transmitted the ranging signal, and (iv) an updating mechanism for updating the initial 
location estimate, using the at least one parameter that is related to the range to a 
beacon of one of the pluralities of beacons, thereby providing a first updated estimate 
of the location of the mobile unit, and for updating the first updated estimate, using 
the at least one parameter that is related to the range to a beacon of another of the 
1 5 pluralities of beacons. 

In the broadest sense, the present invention includes methods and systems that 
enable a mobile unit to determine its location by appropriate processing of signals 
received from mobile beacons and stationary beacons. As understood herein, a 
beacon is a device that includes a transmitter that transmits a signal, typically a radio 
20 frequency signal, that the mobile unit receives and then uses to determine its own 
location. Some of the beacons are stationary, meaning that these beacons are in fixed 
and known locations. Other beacons are mobile, meaning that these beacons move 
along known trajectories, such that at any time, the instantaneous locations of the 
relevant mobile beacons are known. Typically, the mobile unit is a cellular telephone, 
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the stationary beacons are base transceiver stations of a cellular telephone network, 
and the mobile beacons are satellites of a low-Earth-orbit-satellite-based positioning 
system such as the GPS system and the Russian GLONASS system, the trajectories of 
the mobile beacons then being the orbits of the satellites, as described by the satellite 
5 ephemeris and time models. More generally, however, the scope of the present 
invention includes any suitable beacons. For example, the stationary beacons may be 
special purpose beacons set up specifically to enable the tracking of a mobile unit, and 
the mobile beacons may be aircraft flying along known trajectories. 

Although the scope of the present invention extends to the use of any suitable 
10 types of mobile and stationary beacons, the present invention is described herein in 
terms of the preferred and most typical mobile and stationary beacons: GPS satellites 
and base stations (BTSs) of a cellular telephone network. It will be clear to those 
skilled in the art how to adapt the principles of the present invention to other types of 
mobile and stationary beacons. 
15 The mobile unit receives timing signals from the mobile and stationary 

beacons. In the case of the mobile beacons being GPS satellites, the timing signals 
are successive PN frames, modulated by the Satellite Data Messages. In the case of 
the stationary beacons being BTSs, the timing signals are synchronization burst 
sequences transmitted in the cellular network control channel. 
20 In order for these signals to be interpreted correctly by the mobile unit, if the 

cellular telephone network is asynchronous, with the BTSs being synchronized neither 
to the GPS network nor to each other, then the time offsets between the individual 
BTSs and the GPS system must be known by the mobile unit. Therefore, at least one 
reference unit is provided, for determining offsets between the BTSs and GPS time. 
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One reference unit is associated with one or more BTSs. Each reference unit includes 
a GPS-based time reference module for receiving the PN frames from the GPS 
satellites and cross-correlating them with local copies of the PN sequences to 
determine GPS time, a cellular communications module for receiving the 
5 synchronization burst sequences from the associated BTSs, and a synchronization 
module for determining the timing offsets of the associated BTSs relative to GPS 
time. The cellular telephone network is provided with a locator service center, and the 
reference units transmit the timing offsets to the mobile unit via the locator service 
center. Optionally, one of the reference units is designated as a master reference unit 
10 that also transmits the Satellite Data Messages to the mobile unit via the locator 
service center. (This master reference unit may be physically integrated with the 
locator service center.) The mobile unit is synchronized to the cellular network in the 
conventional manner. This allows the mobile unit to also be synchronized 
approximately to the GPS network: having received the BTS timing offsets from the 
15 reference units via the locator service center, the mobile unit computes the time offset 
between its own clock and GPS time to within about 30 microseconds. 

The mobile unit receives the satellite timing signals collectively from the GPS 
satellites that are in view, digitizes the combined signal, and applies matched filter 
algorithms to the combined signal to extract pseudoranges and rates of change of 
20 pseudoranges from the combined signal, as follows. Successive frames of the 
digitized combined signal are arranged as columns of a matrix. As many instances of 
the matrix are formed as there are satellites in view. That the mobile unit is 
approximately synchronized with the GPS network allows each matrix instance to be 
demodulated, with respect to the Satellite Data Message of the corresponding satellite, 
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by multiplying groups of columns of the matrix instance by corresponding bits of the 
Satellite Data Message of the corresponding satellites. If a master reference unit is 
provided, then any PN frames may be used to form the instances of the matrix; 
otherwise, only those PN frames, that correspond to bits of the Satellite Data 

5 Messages that are common to all satellites, are used. The rows of the matrices are 
transformed using a discrete orthogonal transform, preferably a discrete Fourier 
transform. The elements of the matrices are multiplied by respective Doppler 
compensation factors. The columns of the matrices are convolved with the PN 
sequences of the in-view satellites to provide, for each in-view satellite, a signal 

10 ambiguity function matrix. The squared magnitudes of the elements of these matrices 
include local peaks, the row and column coordinates of which correspond to candidate 
hypotheses for the pseudoranges to the corresponding satellites and to the rates of 
change of these pseudoranges. This procedure is repeated several times to provide 
several signal ambiguity function matrices for each in-view satellite. Then, for each 

15 in-view satellite, the best hypothesis is chosen, for example, by integrating the 
matrices noncoherently and then identifying the most likely peak in the sum; or by 
passing the coordinates of all the peaks of the matrices to a "track before detect" 
algorithm. In either case, the peaks are identified using a constant false alarm rate 
detection scheme. 

20 The pseudoranges to stationary beacons are simply the differences between the 

measured times of arrival of the synchronization burst sequences and the transmission 
times of the synchronization burst sequences according to the GPS system clock, 
multiplied by the speed of light. Note that, because the mobile unit has received the 
time offsets between the BTS clocks and the GPS system clock, from the reference 


units via the locator service center, the mobile unit can refer these transmission times 
to the GPS system clock. Note also that, as a consequence, the pseudoranges to the 
BTSs are with reference to the same range offset, c b , as applies to the GPS satellites. 
As each stationary beacon pseudorange is acquired by the mobile unit, and as each 
mobile beacon pseudorange - pseudorange rate pair is acquired by the mobile unit, 
these data are input as update measurements to a dynamic estimator to produce an 
updated estimate of the location of the mobile unit. Preferably, the dynamic estimator 
is based on a linear optimal estimator, such as a Kalman filter. The dynamic estimator 
requires an initial estimate of the location of the mobile unit. This initial estimate is 
obtained by triangulation. Alternatively, under certain conditions, the location of a 
predefined reference point, such as the nearest stationary beacon, may be used as the 
initial estimate of the mobile unit's location. 

BRIEF DESCRIPTION OF THE DRAWINGS 

The invention is herein described, by way of example only, with reference to 
the accompanying drawings, wherein: 

FIG. 1 is a block diagram of a mobile unit of the present invention; 

FIG. 2 is a schematic diagram of a locator system of the present invention; 

FIG. 3 is a block diagram of a reference unit of the present invention; 

FIG. 4 is a flow chart of the joint estimation of position from satellite 
pseudoranges, rates of change of satellite pseudoranges and base station pseudoranges. 
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DESCRIPTION OF THE PREFERRED EMBODIMENTS 

The present invention is of a system and method for locating a mobile unit, for 
example a mobile unit of a wireless communications system, based on timing signals, 
both from mobile beacons such as GPS satellites and from stationary beacons such as 
5 base stations of the wireless communications system. Specifically, the present 
invention can be used to locate a suitably equipped cellular telephone in an 
environment, such as an urban environment, in which it is difficult to establish 
simultaneous lines of sight to four or more GPS satellites. 

Although the present invention is described herein in reference to the GPS 
10 system and a cellular telephony network, it will be appreciated that this description is 
purely illustrative, and that the scope of the present invention extends to the use of any 
suitable types of mobile and stationary beacons, as defined above. 

The principles and operation of a locator system and method according to the 
present invention may be better understood with reference to the drawings and the 
15 accompanying description. 

Referring now to the drawings, Figure 1 is a partial high level block diagram 
of a mobile unit 10 of the present invention. Mobile unit 10 includes an antenna 12 
for receiving RF signals transmitted by GPS satellites, a downconverter 14 for 
downconverting the RF signals received from the GPS satellites to an intermediate 
20 frequency (IF), a digitizer 16 for digitizing the IF signals, a clock 22 for providing 
local timing signals, a memory 18 for storing the digitized GPS signals , and a GPS 
processor 20 for processing the signals stored in memory 18. Mobile unit 10 also 
includes an antenna 24 for exchanging RF signals with BTSs of a cellular telephone 
network, a cellular modem 26 for downconverting and digitizing the control channel 
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RF signals received from the BTSs, a memory 27 for storing the digitized BTS 
signals, and a cellular processor 28 for processing the digitized BTS signals with 
reference to timing signals from clock 22, as described below. The outputs of GPS 
processor 20 and cellular processor 28 are reconciled in a joint estimator 30, as 
5 described below, to produce a joint estimate of the location of the mobile unit of 
which mobile unit 10 is a part. GPS processor 20, cellular processor 28 and joint 
estimator 30 function under the overall supervision of a management unit 29. In 
particular, management unit 29 formulates requests for supplemental data, and, using 
cellular modem 26 and antenna 24, transmits these requests as RF signals to the BTSs, 

10 as described below. 

Figure 2 is a schematic diagram of a locator system 40 of the present 
invention. Locator system 40 is based on a conventional cellular telephone network 
42 that includes multiple BTSs 44 , and on the GPS system, represented in Figure 2 by 
satellites 70. The terrestrial component of locator system 40 adds to cellular 

15 telephone network 42 a locator service center 46 and a multitude of reference units 50, 
at fixed and known locations. Each reference unit 50 is associated with one or more 
BTSs 44, and communicates with those BTSs 44 and, via the nearest of those BTSs 
44, with service center 46. Mobile unit 10, using antenna 24 and cellular modem 26, 
also receives messages from service center 46 via the nearest BTS 44, as described 

20 below. 

Figure 3 is a high level block diagram of reference unit 50. Reference unit 50 
includes a cellular communications module (CCM) 52 for communicating with 
service center 46 via cellular network 42, a GPS-based time reference module 66, and 
a synchronization module 60 for synchronizing CCM 52 to GPS-based time reference 
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module 66 as described below. CCM 52 includes both transmission functionality 54 
and reception functionality 56, via an antenna 58. GPS-based time reference module 
66 receives GPS satellite signals via an antenna 68. Synchronization module 60 
includes both a clock 64 for providing reference timing signals and a processor 62 for 
5 establishing timing offsets between the GPS system and the BTSs 44 that are 
associated with reference unit 50, as described below. The main function of each 
reference unit 50 is to determine these timing offsets and to transmit these timing 
offset to locator service center 46. Optionally, one of reference units 50 is designated 
as a master reference unit. The master reference unit has the additional function of 
10 receiving the Satellite Data Messages from satellites 70 and transmitting these 
Satellite Data Messages to locator service center 46. 

Locator service center 46 is a node of cellular telephone network 42 whose 
function is to coordinate communications between mobile unit 10 and cellular 
telephone network 42 that relate to the self-location functionality of mobile unit 10. 
15 In particular, locator service center 46, upon receipt of a query from mobile unit 10, 
transmits to mobile unit 10 supplemental data such as the timing offsets between the 
GPS system and one or more BTSs 44 and optionally the Satellite Data Messages of 
in-view satellites 70. Optionally, the functionality of joint estimator 30 may reside in 
location service center 46, with the outputs of GPS processor 20 and cellular 
20 processor 28 being transmitted to location service center 46 for reconciliation. 

The signal received by mobile unit 10 from the GPS system, as a function of 
time t, is a sum of the signals G(t) transmitted by GPS satellites 70 that are in view at 
time The signal received by mobile unit 10 from any one satellite 70, indexed by an 
index j, is in the form 
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Gj{i) = KjDjit-x^gjit-x^expliGyit-Xj)] 
where / is the square root of -1 ; t, is the true one-way propagation time of radio waves 
from the satellite to mobile unit 10; K, is an amplitude factor that depends on the true 
range Ry=cXj to satellite j at time t, and on the GPS antenna gain pattern; Dj{t-x } ) is the 
Satellite Data Message of satellite j at time t-xy gj(t-Xj) is the PN sequence of satellite 7 
at time t-xy and exp[ia(t-Xj)] is a Doppler factor whose Doppler frequency <o depends 
on the radial component of the velocity of satellite/ relative to mobile unit 10, on the 
bias of clock 22, and on other imperfections in mobile unit 10. The total signal G(f) 
received by mobile unit 10 is the sum of these signals, ^Gj(t). 

The object of the processing in GPS processor 20 is to measure py=Rj+c b and 
its rate of change dpjdt for all the satellites in view. Note that all the pseudoranges p, 
are referred to the same range offset c h . As in the prior art methods, c b is one of the 
unknowns that is solved for subsequently to determine the location of mobile unit 10, 
so it is not necessary to measure x, explicitly. 

As in US 5,663,734, several hundred seconds worth of a digitized total GPS 
signal is stored in memory 18 and then processed by GPS processor 20. One 
difference between US 5,663,734 and the GPS processing of the present invention is 
that according to the present invention, groups of many more than 5 to 10 PN 
sequence frames are processed coherently. 

The total signal received by mobile unit 10 is oversampled by digitizer 16, 
starting at an initial reception time t 0 , typically at a sampling rate of between two and 
four samples per chip. With four samples per chip, there are a total of 4x1023=4092 
samples per frame. For each in-view satellite 70, indexed by index;', these samples 
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are stored as successive columns of a matrix My that, in the case of a sampling rate of 
four samples per chip, has 4092 rows and as many columns as the number of 
consecutive frames included in the computation: the first column of My, matrix 
elements M^U) through A^(4092,l) are the first 4092 samples, corresponding to the 
5 first millisecond's worth of data; the second column of My, matrix elements M^l,!) 
through ^(4092,2), are the second 4092 samples, corresponding to the second 
millisecond's worth of data; etc., through a number of columns n equal to the number 
of frames of data that are stored in memory 18, at one millisecond per frame. 
Preferably, between 100 milliseconds and 1000 milliseconds worth of data (100 to 
10 1000 frames) are stored in memory 18 as matrix My for each in-view satellite 70. 

The first processing step is to demodulate each matrix My, i.e., to remove the 
influence of Satellite Data Message D } . This is done by multiplying groups of 
columns of My by the corresponding bits of D,(t-Xj). A necessary prerequisite for this 
demodulation is for the bits of D } to be aligned with the columns of My to within 
15 about one millisecond. This alignment is referred to herein as "approximate 
synchronization" of mobile unit 10 with the GPS system. 

Cellular telephone network 42 may be synchronized to the GPS system clock 
in any case, for example, if cellular telephone network operates according to the 
CDMA IS-95 standard. Mobile unit 10 is synchronized to cellular telephone network 
20 42 in the conventional manner. Because cellular telephone network 42 is time-aligned 
(synchronized) with GPS system time, mobile unit 10 also is synchronized to GPS 
system time and therefore is able to demodulate matrix My relative to Satellite Data 
Message D,. 
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If cellular telephone network 42 is not synchronized to the GPS system clock, 
and, in particular, if BTSs 44 are not synchronized to each other, mobile unit 10 
nevertheless can obtain, from reference units 50 via locator service center 46, the 
transmission times of the synchronization burst sequences of BTSs 44 according to 
5 the GPS clock. When CCM 52 of a reference unit 50 receives a synchronization burst 
sequence from one of the BTSs 44 associated with reference unit 50, synchronization 
module 60 determines the time offset between that BTS 44 and the GPS system with 
reference to clock 64. Because the locations of reference unit 50 and the associated 
BTSs 44 are fixed and known, the propagation time from each associated BTS 44 to 

10 reference unit 50 also is known, so the time offset is referred to the time at which that 
BTS 44 transmitted the synchronization burst sequence, and not the time at which 
reference unit 50 received the synchronization burst sequence. CCM 52 then 
transmits this time offset to locator service center 46. This is done for the 
synchronization burst sequence of each associated BTS 44, so that locator service 

15 center 46 always has a current list of time offsets of all BTSs 44 from the GPS system. 
Finally, mobile unit 10 receives the current time offsets of the relevant BTSs 44 from 
locator service center 46 on request. The time offset of the nearest BTS 44 from the 
GPS system is used by mobile unit 10 to correct the nominal transmission times of the 
synchronization burst sequences from the nearest BTS 44 for this offset, and so to 

20 obtain an estimate of its own time offset To from GPS system time to within 30 
microseconds. The time offsets of BTSs 44 from the GPS system also are used to 
obtain BTS pseudoranges that are referred to range offset c/,, as described below. 

If the frames that are used to construct matrix Mjj correspond to Satellite Data 
Message header bits or to Satellite Data Message handover word bits, which are 
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common to the Satellite Data Messages of all satellites 70, then mobile unit 10 either 
knows these bits a priori (in the case of header bits) or computes these bits based on 
the time of the week (in the case of handover word bits). In the embodiment of the 
present invention in which a master reference unit is designated to receive the Satellite 
5 Data Messages and relay these Satellite Data Messages to locator service center 46, 
mobile unit 10 sends a request to locator service center 46 to obtain Satellite Data 
Message D, so that the satellite-dependent bits of Satellite Data Message Dj can be 
used to demodulate matrix Mjj. 

The second processing step is to transform each matrix Afj, into a matrix M B j 
10 by performing a discrete orthogonal transform on each of the rows of matrix My, 
thereby transforming the rows from the time domain to the frequency domain, so that 
matrices Mnj are in a mixed time-frequency domain. The most preferred discrete 
orthogonal transform is the discrete Fourier transform, although the scope of the 
present invention includes other transforms, such as the Walsh-Hadamard transform, 
1 5 the various kinds of Wavelet transforms and the Hartley transform. 

The third processing step is the multiplication of the elements of each matrix 
M rij by Doppler compensation factors. The form of the Doppler compensation factor 

for the (r,d) element of matrix M nj isex p( _ ~f"("~ + ^) ' where m is the number of 

rows in matrix M Uj , n is the number of columns in matrix M U j and h is the integer part 
20 of the Doppler frequency shift of the signal G/t), measured in KHz. Note that, with 
one one-millisecond PN frame per column of matrix My, each row of matrix M uj 
spans one KHz. Initially, the true Doppler frequency shift of each satellite 70 can be 
estimated by mobile unit 10 to within about 4 KHz. Therefore, four different sets of 
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Doppler compensation factors are needed, resulting in four different matrices M U j for 
each satellite j. Subsequently, as described below, much more accurate estimates of 
the Doppler frequency shifts of satellites 70 are obtained, and only one matrix M U j is 
needed per satellite For clarity, the remaining discussion of how GPS pseudoranges 
and pseudorange rates of change are measured according to the present invention is in 
terms of a single matrix M U j per satellite/ 

The fourth processing step is to convolve each of the columns of each matrix 
Muj with the PN sequence of satellite j, to produce a matrix Mj that is equivalent to 
the matched filter, in range and frequency, known in the field of radar signal 
processing as the "signal ambiguity function". See, for example, Fred E. Nathanson, 
Radar Design Principles, Second Edition (McGraw-Hill, 1991). It is this use of 
demodulated signal ambiguity functions, subsequent to approximate synchronization 
of mobile unit 10 to the GPS system, that enables the present invention to process 
coherently many more than the 5 to 10 PN sequence frames that are processed 
coherently by Krasner in US 5,663,734. The preferred convolution method is the 
well-known fast convolution method, also used by Krasner in US 5,663,734: perform 
a discrete Fourier transform on a column of matrix Mnj, multiply the transformed 
column by the discrete Fourier transform of the PN sequence, and invert the discrete 
Fourier transform of the column. Note that the multiplication of the elements of 
matrix M m by the satellite-dependent parts of the Doppler compensation factors, 


transformed column of M nj , appropriately shifting the discrete Fourier transform of 
the PN sequence prior to the multiplication of the transformed column by the discrete 



, may be deferred to this processing step and implemented by, for each 
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Fourier transform of the PN sequence. Note also that the matrix elements of Mj are 
complex. The overall appearance of matrix M h in terms of the squared magnitudes of 
its elements, is a set of peaks above a general background level. The row coordinate 
of a peak is its pseudorange (units of distance) and the column coordinate (units of 
5 frequency) of a peak is proportional to the rate of change of its pseudorange (units of 
velocity). 

Although the preferred structure of matrix My has one frame per column, it 
also is possible to construct matrix M V] with more than one frame per column, as long 
as an integral number of frames are used per column. Using one frame per column 

10 gives the optimal resolution for determining the rates of change of the pseudoranges. 

A constant false alarm rate detection scheme is used to identify the peaks. 
Preferably, the detection scheme is based on an adaptive threshold, based on a 
background level that is obtained, matrix element by matrix element, by averaging the 
squared magnitudes of all the matrix elements in a window (for example, an 8x8 

15 window or a 16x16 window) around the target matrix element. The set of matrix 
elements, whose squared magnitudes exceed their local average by a threshold factor, 
is in the form of isolated clusters of matrix elements, with each cluster corresponding 
to one peak. The coordinates of the peaks are the centroids of the clusters, weighted 
by the magnitudes of the matrix elements in the clusters. The threshold factor is 

20 determined by the required maximum false alarm rate, as is well-known in the art. 

The output of this procedure is a set of hypotheses for the pseudorange p, and 
its rate of change dpjdt associated with the initial reception time f 0 . At this point, the 
hypothesis associated with the peak most likely to be the most reliable measure of pj 
and dpjdt, for example, the highest peak, could be selected to provide the coordinates 
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that determine p 3 and dpj/dt. Preferably, however, several instances of the matrices 
Mj, based on several successive instances of the digitized total GPS signal, are 
generated. Preferably, between five and ten such instances of each matrix Mj are so 
generated. Subsequent processing is directed towards obtaining, for satellite j, values 

5 of py and dpjdt that are based on all the instances of Mj processed collectively. 

One method of processing the instances of Mj collectively is to integrate them 
noncoherently. In standard noncoherent integration, the matrix elements of all the 
instances are replaced by their squared magnitudes; the instances are added to each 
other, matrix element by matrix element, to produce a sum matrix; the false alarm 

10 detection scheme described above is used to identify the peaks of the sum matrix; and 
the coordinates of the peak most likely to be the most reliable measure of pj and 
dpjldt, for example, the highest peak, are taken as the values of p, and dpjdt. In 
binary noncoherent integration, the matrix elements of all the instances are replaced 
by their squared magnitudes, but then the constant false alarm detection scheme is 

1 5 used to compare each matrix element with a threshold. Matrix elements that exceed 
the threshold are replaced by ones, and all other matrix elements are replaced by 
zeros. The instances now are added to each other, matrix element by matrix element, 
to produce a sum matrix, and a second threshold is applied to the sum matrix to 
identify its peaks. Again, the coordinates of the peak, of the sum matrix, that is most 

20 likely to represent the most reliable measure of p 7 and dpjdt, are used as the values of 
Pj and dpjldt. 

Another way to obtain values of p 7 and dpjdt that are based on all the 
instances of Mj processed collectively is to replace the matrix elements of all the 
instances with their squared magnitudes, as in noncoherent integration, but then to 
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apply the constant false alarm rate detection scheme separately to each instance, to 
obtain, for each instance, a set of hypotheses for p, and dpjdt. These hypotheses and 
their associated start times are input to a "track before detect" identification scheme to 
identify the most likely hypothesis. "Track before detect" is a method that is well 
5 known in radar technology for identifying targets in a low signal-to-noise 
environment. See, for example, the book by Fred E. Nathanson cited above. 

The pseudoranges to BTSs 44 are the differences between the measured times 
of arrival of the synchronization burst sequences and the nominal transmission times 
of the synchronization burst sequences according to the GPS system clock. If cellular 
10 telephone network 42 is synchronized to the GPS system clock, then the nominal 
transmission times of the synchronization burst sequences, according to clock 22, are 
offset from their true times by the same time offset, 7b, as the time offset by which the 
nominal transmission times of the PN sequences by satellites 70, according to clock 
22, are offset from their true times. As a result, all the pseudoranges, both to satellites 
15 70 and to BTSs 44, share the same range offset c b . If cellular telephone network 42 is 
not synchronized to the GPS system clock, then mobile unit 10 obtains the time 
offsets of BTSs 44 from the GPS system from locator service center 46, as described 
above in the context of the approximate synchronization of mobile unit 10 with the 
GPS system for the purpose of demodulation. Mobile unit 10 then uses these time 
20 offsets to correct the times of arrival of the synchronization burst sequences from 
BTSs 44, as measured by clock 22, so that the pseudoranges to BTSs 44 are referred 
to the same range offset c b as the pseudoranges to satellites 70. The pseudoranges to 
the in-view satellites 70, the rates of change of the pseudoranges to the in-view 
satellites 70, and the pseudoranges to BTSs 44 are input to joint estimator 30. 
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The method used by joint estimator 30 to combine pseudoranges p and 
pseudorange rates of change dp/dt from GPS processor 20, and estimated 
pseudoranges p from cellular processor 28 to derive the location of mobile unit 10 
now will be discussed. The essence of this method is to use each new measurement (a 
5 pseudorange p and pseudorange rate of change dp/dt relative to a particular satellite 
70 or a pseudorange p relative to a particular BTS 44), as that measurement is 
obtained, to update an a priori estimate of the location vector r=(x,y^) T and of the 
velocity vector v={v x ,v y ,v z ) T of mobile unit 10. (It is for this reason that this 
measurement is referred to herein as an "update measurement".) If mobile unit 10 
10 were to remain stationary between successive measurements and if range offset c b 
were to remain constant between successive measurements, or if four or more 
measurements were to be obtained from matrices Mj computed from signals received 
simultaneously by mobile unit 10 from the GPS system, it would be possible to fix the 
location of mobile unit 10 by triangulation. Because mobile unit 10 moves and c b may 
15 drift between successive measurements, if measurements obtained from GPS signals 
received by mobile unit 10 at different times are to be used, or if pseudoranges to 
BTSs 44 are to be used, a method of estimating r that accounts for the motion of 
mobile unit 10 and for the drift of clock 22 must be used. Such a method is referred to 
herein as a "dynamic" estimation method. The preferred dynamic estimation methods 
20 belong to the class of linear optimal estimators. The most preferred dynamic 
estimation method is based on a Kalman filter. This method now will be described in 
detail. 

In-between updates of r and v, mobile unit 10 is assumed to move in a known 
manner. The model of the motion of mobile unit 10 used herein as an example is 
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linear motion with additive random noise. Other models, such as motion constrained 
to be along a known road, also may be used. According to the linear model of the 
motion of mobile unit 10, if the location vector of mobile unit 10 at time t kA is r kA and 
the velocity vector of mobile unit 10 at time 4-i is v k . u then the location and velocity 

5 vectors at time t k =t kA +8 k are r k =r kA +v k .&k and v k =v kA . The range bias c b is defined as 
the difference between the time measured by clock 22 and the GPS clocks, multiplied 
by the speed of light c, and the range drift c d is defined as d(c b )/dt. Then c b(k) =c b(k . 
D+c d(kA) 8 k and c m =c d( k-\). The (eight component) navigation vector n is defined as 
the concatenation of r, v, c b and c d : n=(x,y^,v x ,v y ,v z ,c b ,c d ) T . 

10 The state vector eF=i&c,8yJ5z,8v x ,bvy,8v zt 8c b ,8c d ) T of the Kalman filter is the 

estimated navigation error, i.e., the difference between the true navigation vector and 
the linearly propagated navigation vector. This state vector obeys the propagation 
equation 

e k =A k -\e k -i + w k -i 
1 5 where the transition matrix A k .\ is: 
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and w kA is the process noise. w k is defined via the model noise covariance matrix 
Gm: 


Qk-\ = w kA w kA T = 


'RW X 0 0 0 0 0 0 0" 

0 RW y 0 0 0 0 0 0 

0 0 000000 

0 0 000000 

0 0 000000 

0 0 000000 

0 0 000000 
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where RW X and RW y are the covariances of the process noise in the x and y directions 
(note that mobile unit 10 is assumed to move substantially in the (x,y) plane), modeled 
as random walks, and Wis a range drift error covariance, modeled as white noise. 

Immediately following an update of the location vector r and of the velocity 
vector v of mobile unit 10, the state vector e is zero by definition. Subsequently, until 
the next measurement, as the navigation vector n is propagated forward in time, the 
covariance of the state vector, an 8 x 8 matrix P, also is propagated forward in time: 
P k =AP + kA A T +Q6 k 

Here, the superscript "+" refers to the a postiori estimate of P after the most recent 
update, and the superscript "-" refers to the a priori estimate of P, before the next 
update. The initialization of P is discussed below. 

A measurement obtained from a signal from a satellite 70 is a two component 
vector, z={p,dpldtf. The matrix H h that relates the state vector e to this measurement 
vector at time t k is: 


where the vector u=(u x ,u y ,u z ) is a unit vector from the estimated location of mobile 
unit 10 to satellite 70: 



\ u y u z 0 0 0 1 0 s 
0 0 0 u x u y u z 0 1, 


u = (r k -s k )l\rk-Sk\ 


r k is the a priori mobile unit 10 location vector, obtained by forward propagation as 
described above. s k is the position vector of satellite 70 at time t k , obtained by GPS 
processor 20 using the satellite ephemeris and time model of satellite 70. 

A measurement obtained from a signal from a BTS 44 is a scalar (a one- 
5 component vector), z=p. The matrix H k that relates the state vector e to this 
measurement vector at time t k is: 

H k = (u x u y u z 0 0 0 1 0) 
where the vector u=(u x ,u y ,u z ) is a unit vector from mobile unit 10 to BTS 44. Note 
that the positions of BTSs 44 are fixed and known. 
10 In both cases, an estimate of e is obtained by multiplying the measurement by 

the Kalman gain matrix K k : 

K k = F k H k T {H k F k H k T +R k ) A 
e = K k z 

where the measurement noise covariance matrix R k is: 



for a satellite measurement and 

for a BTS measurement. The diagonal elements of a k are the variances of the 
corresponding pseudorange and pseudorange rate of change measurement errors. A 
20 new a postiori estimate of P is obtained as follows: 
P\ = (I-K k H k )F k 
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where / is the unit matrix. Finally, the a postiori navigation vector n k is obtained by 
adding the state vector e to the a priori navigation vector n k : x k + =x k +bx, y k =y k +by, 
z k =z k +bz, v xk =v xk +hv x , v yk =v yk +hv y , v zk + =v zk +8v z , c bk + =c bk +5c b and c dk + =c dk +8c d . 

Figure 4 is a flow chart of the overall processing effected by joint estimator 30 
5 to track the motion of mobile unit 10. The processing is initialized (block 100) with 
an initial estimate of the navigation vector n, and by setting the diagonal elements of 
the covariance matrix P equal to initial values that reflect the initial standard 
deviations of the errors in r, v, c b and c d . How the initial estimate of n is obtained is 
discussed below. In block 102, joint estimator 30 determines whether a measurement 
10 vector z has been obtained from GPS processor 20 or from cellular processor 28. If a 
measurement has not been obtained, then joint estimator 30 propagates (block 104) 
the navigation vector n forward in time, as described above. If a measurement has 
been obtained, then joint estimator uses the Kalman filter described above (block 106) 
to obtain an a postiori estimate of the navigation vector n. Note that this updating of 
15 n is based on a single new measurement, in contrast to the use of a Kalman filter in 
prior art GPS navigation, in which at least four independent simultaneous 
measurements are needed. This a postiori estimate is reported (block 108), and joint 
estimator 30 again checks for a new measurement vector z from GPS processor 20 or 
from cellular processor 28. 
20 There are several methods by which the initial estimate of the navigation 

vector n may be obtained. Most simply, a static estimate is used as the initial estimate 
of r in block 100, and v is initially set to zero. By "static" estimation is meant that the 
estimation is performed under the assumption that mobile unit 10 did not move while 
the measurements were made and that the range offset c b also was stationary while the 
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measurements were made. The simplest "static" estimate of r is the location of a 
predefined reference point, such as the BTS 44 nearest to mobile unit 10. c b and c d 
also are initialized to zero. This estimate is suitable for use in block 100 of Figure 4 if 
the first four measurements are satellite measurements, because the error in using the 
5 location of a terrestrial reference point is much smaller than the distances to the 
satellites. Triangulation is used to combine initial measurements from satellites with 
initial measurements from BTSs 44, with the proviso that the equation used with a 
measurement from a BTSs 44 is the full (quadratic) triangulation equation 
\s-r\ = p-c b 

10 rather than the linearized equation because BTSs 44 are too close to mobile unit 10 to 
justify linearization of this equation. Mathematical techniques for solving 
simultaneous quadratic and linear equations are well-known and need not be 
elaborated here. 

Once the tracking of mobile unit 10 has commenced, the processing described 
15 above may be made both more accurate and more efficient. For example, the Doppler 
frequency shifts of the signals G/f) may be estimated more accurately. The total 
Doppler frequency shift is the sum of three terms: the Doppler frequency shift due to 
the motion of satellite J, the Doppler frequency shift due to the motion of mobile unit 
10, and a contribution from the drift of clock 22. The first term can be computed from 
20 the ephemeris and time model of satellite j. Mobile unit 10 moves sufficiently slowly 
that the second term is negligible. The main initial uncertainty in the total Doppler 
frequency shift lies in the third term; but this term may be inferred from the range drift 
component c d of the navigation vector «. It then is possible to obtain an estimate of h, 
the integer part of the Doppler frequency shift of G,(t), that is much more accurate 
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than the initial estimate described above, for use in the Doppler compensation factors 
for matrix M Uj . Furthermore, the columns of Mj that are most likely to contain the 
peak corresponding to p y and dp/dt can be predicted to within 10 to 20 Hz, as opposed 
to the 1 KHz range of the full matrix Mj, so that only these columns of matrix Mnj 

5 need to be produced by the discrete orthogonal transform of matrix My. 

Alternatively, having computed a pseudorange and a rate of change thereof for 
one satellite 70, mobile unit estimates the clock drift term of the total Doppler 
frequency shift by subtracting the computed Doppler frequency shift of the nearest 
BTS 44 from the apparent Doppler frequency shift that corresponds to the rate of 

10 change of the pseudorange to this one satellite 70. This estimated clock drift term is 
added to the known Doppler frequency shifts due to the motion of the other satellites 
70 to provide estimates of the total Doppler frequency shifts of the other satellites 70. 

While the invention has been described with respect to a limited number of 
embodiments, it will be appreciated that many variations, modifications and other 

1 5 applications of the invention may be made. 
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WHAT IS CLAIMED IS: 

1 . A method for determining a pseudorange, and a rate of change thereof, 
to a beacon that transmits a signal including a plurality of blocks, each block including 
a plurality of frames of a pseudonoise sequence, each block being multiplied by a bit 
of a data sequence, the signal being shifted relative to a nominal frequency by a 
frequency shift, the method comprising the steps of: 

(a) receiving the signal; 

(b) digitizing said received signal, thereby producing a digitized signal 
including a plurality of bits; 

(c) arranging said digitized signal as columns of an input matrix that 
includes a plurality of rows, each said column including consecutive 
said bits of said digitized signal that correspond to an integral number 
of frames of the pseudonoise sequence; and 

(d) performing a discrete orthogonal transform on each said row of said 
input matrix, thereby producing a transformed matrix. 

2. The method of claim 1 , wherein said integral number is 1 . 

3. The method of claim 1, wherein said discrete orthogonal transform is a 
discrete Fourier transform. 

4. The method of claim 1, wherein said transformed matrix includes a 
plurality of elements, the method further comprising the step of: 
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(e) multiplying each said element of said transformed matrix by a 
respective Doppler compensation factor. 


5. The method of claim 4, further comprising the step of: 
(f) providing an estimate of the frequency shift; 

said Doppler compensation factors being based on said estimate of the frequency shift. 

6. The method of claim 5, wherein said transformed matrix has fewer said 
columns than said input matrix, said estimate of the frequency shift being used to 
determine which said columns of said transformed matrix are produced by said 
discrete orthogonal transform. 

7. The method of claim 4, further comprising the step of: 

(f) convolving each said column of said transformed matrix with the 
pseudonoise sequence. 

8. The method of claim 7, further comprising the step of: 

(g) identifying at least one peak in said transformed matrix, each said at 
least one peak having a row coordinate and a column coordinate. 

9. The method of claim 8, further comprising the step of: 

(h) selecting, from among said at least one peak, a most likely peak, said 
row coordinate of said most likely peak then corresponding to the 
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pseudorange, and said column coordinate of said most likely peak then 
corresponding to the rate of change of the pseudorange. 

10. The method of claim 8, wherein said receiving, said digitizing, said 
arranging, said performing, said multiplying, said convolving and said identifying are 
effected on a plurality of instances of the transmitted signal, the method further 
comprising the step of: 

(h) inferring the pseudorange and the rate of change of the pseudorange 
from said row coordinates and said column coordinates. 

11. The method of claim 10, wherein said inferring is effected using a 
track-before-detect identification scheme. 

12. The method of claim 7, wherein said receiving, said digitizing, said 
arranging, said performing, said multiplying and said convolving are effected on a 
plurality of instances of the transmitted signal to produce a corresponding plurality of 
said transformed matrices, the method further comprising the steps of: 

(h) integrating said transformed matrices non-coherently to produce a sum 
matrix; and 

(i) identifying at least one peak in said sum matrix, each said at least one 
peak having a row coordinate and a column coordinate. 


13. The method of claim 12, further comprising the step of: 

(j) selecting, from among said at least one peak, a most likely peak, said 
row coordinate of said most likely peak then corresponding to the 
pseudorange, and said column coordinate of said most likely peak then 
corresponding to the rate of change of the pseudorange. 

14. The method of claim 4, further comprising the step of: 

(g) providing a plurality of estimates of the frequency shift, thereby 
producing, for each said element of said transformed matrix, a 
corresponding plurality of said respective Doppler compensation 
factors that are based on said estimates of the frequency shift; 
and wherein said multiplying is effected separately for each said estimate of the 
frequency shift, thereby producing a corresponding plurality of Doppler-compensated 
matrices. 

1 5. The method of claim 1 , further comprising the step of: 

(e) multiplying groups of said bits of said digitized signal, that correspond 
to the blocks of the transmitted signal, by respective bits of the data 
sequence, prior to said arranging of said digitized signal as columns of 
said input matrix. 

1 6. A method for determining a pseudorange, and a rate of change thereof, 
to each of a plurality of beacons, each beacon transmitting a respective signal 
including a plurality of blocks, each block including a plurality of frames of a 
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pseudonoise sequence, each block being multiplied by a bit of a data sequence, the 
pseudonoise sequences and the data sequences being uniquely associated with 
respective satellites, all the pseudonoise sequences being of equal length, the 
pseudonoise sequences being mutually orthogonal, the method comprising the steps 
5 of: 

(a) receiving the transmitted signals collectively as a received signal; 

(b) digitizing said received signal, thereby producing a digitized signal 
including a plurality of bits; and 

(c) for each beacon, multiplying groups of said bits of said digitized signal, 
10 that correspond to the blocks of the signal transmitted by said each 

beacon, by respective bits of the data sequence of said each beacon. 

17. The method of claim 16, further comprising the steps of: 

(d) for each beacon, arranging said digitized signal as columns of a 
15 corresponding instance of a first matrix that includes a plurality of 

rows, each said column including consecutive said bits of said digitized 
signal that correspond to an integral number of frames of the 
pseudonoise sequences; and 

(e) for each beacon, performing a discrete orthogonal transform on each 
20 said row of said corresponding instance of said first matrix. 


18. The method of claim 17, wherein each said instance of said first matrix 
includes a plurality of elements, the method further comprising the steps of, for each 
beacon: 
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(f) multiplying each said element of said corresponding instance of said 
first matrix by a respective Doppler compensation factor; and 

(g) convolving each said column of said corresponding instance of said 
first matrix with said pseudonoise sequence of said each beacon. 

19. The method of claim 1 8, further comprising the step of: 

(h) for each beacon, identifying at least one peak in said corresponding 
instance of said first matrix, each said at least one peak having a row 
coordinate and a column coordinate. 

20. The method of claim 1 9, further comprising the step of: 

(i) for each beacon, selecting, from among said at least one peak, a most 
likely peak, said row coordinate of said most likely peak then 
corresponding to the pseudorange of said each beacon, and said column 
coordinate of said most likely peak then corresponding to the rate of 
change of the pseudorange of said each beacon. 

21. The method of claim 20, wherein said receiving, said digitizing, said 
multiplying by respective bits of the data sequence, said arranging, said performing, 
said multiplying of said elements by said Doppler compensation factors, said 
convolving and said identifying are effected on a plurality of instances of the 
transmitted signal, the method further comprising the step of: 

(j) for each beacon, inferring the pseudorange of said each beacon and the 
rate of change of the pseudorange of said each beacon from said row 


coordinates and said column coordinates corresponding to said each 
beacon. 


22. The method of claim 18, wherein said receiving, said digitizing, said 
5 multiplying by respective bits of the data sequence, said arranging, said performing, 
said multiplying of said elements by said Doppler compensation factors and said 
convolving are effected on a plurality of instances of the transmitted signal to produce, 
for each beacon, a corresponding plurality of said second matrices, the method further 
comprising the steps of: for each beacon: 
10 (h) integrating said second matrices non-coherently to produce a sum 

matrix; and 

(i) identifying at least one peak in said sum matrix, each said at least one 
peak having a row coordinate and a column coordinate. 


23. The method of claim 22, further comprising the step of: for each 

beacon 

(j) selecting, from among said at least one peak, a most likely peak, said 
row coordinate of said most likely peak then corresponding to the 
pseudorange, and said column coordinate of said most likely peak then 
corresponding to the rate of change of the pseudorange. 


24. The method of claim 16, wherein said respective bits of said data 
sequences are common to all the beacons. 


25. A receiver for receiving a signal transmitted by a beacon, the signal 
including a plurality of frames of a pseudonoise sequence, comprising: 

(a) an antenna for receiving the transmitted signal; 

(b) a mechanism for digitizing the received signal to produce a digitized 
signal including a plurality of bits; 

(c) a memory for storing the digitized signal as columns of a matrix that 
includes a plurality of rows, each said column including consecutive 
said bits of said digitized signal that correspond to an integral number 
of frames of the pseudonoise sequence; and 

(d) a processor for performing a discrete orthogonal transform on each of 
said rows of said matrix. 

26. A locator system for locating a mobile unit, comprising: 

(a) at least one beacon having a respective pseudonoise sequence and a 
respective data sequence, each said data sequence including a plurality 
of bits, each said at least one beacon operative to transmit a respective 
transmitted signal, each said respective transmitted signal including a 
plurality of blocks, each said block including a plurality of frames of 
said respective pseudonoise sequence, each said block being multiplied 
by a bit of said respective data sequence, 

(b) a reference unit including: 

(i) a reference unit receiver for: 

(A) receiving said at least one transmitted signal as a 
reference unit received signal, and 


(B) recovering said at least one data sequence from said 
respective reference unit received signal, and 
(ii) a transmitter for transmitting said at least one data sequence to 
the at least one mobile unit; and 
5 (c) in the mobile unit: 

(i) a first mobile unit receiver for receiving said at least one data 
sequence, and 

(ii) a second mobile unit receiver including: 

(A) an antenna for receiving the at least one transmitted 
10 signal collectively as a mobile unit received signal, 

(B) a mechanism for digitizing the mobile unit received 
signal to produce a digitized signal including a plurality 
of bits, 

(C) a memory for storing, for each of the at least one 
15 beacon, an instance of said digitized signal, and 

(D) a processor for multiplying each of said at least one 
instance of said digitized signal by said bits of said 
respective data sequence. 


27. The locator system of claim 26, wherein each said at least one beacon 
is operative to move on a known respective trajectory while transmitting said 
respective transmitted signal. 
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28. The locator system of claim 26, wherein said processor is operative, for 
each of the at least one beacon, to arrange said corresponding instance of said 
digitized signal in said memory as columns of a matrix that includes a plurality of 
rows, each said column including consecutive said bits of said corresponding instance 
5 that correspond to an integral number of frames of said respective pseudonoise 
sequence, and to perform a discrete orthogonal transform on each of said rows of said 
matrix. 


29. A method for determining a location of a receiver, comprising the steps 

10 of: 

(a) providing a plurality of beacons having respective pseudonoise 
sequences, all said pseudonoise sequences being of equal length; 

(b) transmitting, by each said beacon, a respective signal including a 
plurality of frames of said respective pseudonoise sequence; 

15 ( c ) receiving said transmitted signals collectively as a received signal, by 

the receiver; 

(d) inferring, for each beacon, a pseudorange and a rate of change of said 
pseudorange; and 

(e) inferring the location of the receiver from said pseudoranges and from 
20 said rates of change of said pseudoranges. 


30. The method of claim 29, wherein said inferring of said pseudoranges 
and said rates of change is effected by steps including: 
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(i) digitizing said received signal, thereby producing a digitized signal 
including a plurality of bits; 

(ii) for each said beacon, arranging said digitized signal as columns of a 
matrix that includes a plurality of rows, each said column including 
consecutive said bits of said digitized signal that correspond to an 
integral number of said frames of said pseudonoise sequences; and 

(iii) for each said beacon, performing a discrete orthogonal transform on 
each said row of said matrix. 

31. A method for determining a pseudorange, and a rate of change thereof, 
to a beacon that transmits a signal including a plurality of blocks, each block including 
a plurality of frames of a pseudonoise sequence, each block being multiplied by a bit 
of a data sequence, the signal being shifted relative to a nominal frequency by a 
Doppler frequency shift, the method comprising the steps of: 

(a) receiving the signal; 

(b) digitizing said received signal, thereby producing a digitized signal 
including a plurality of bits; and 

(c) applying a matched filter algorithm to said digitized signal to extract 
therefrom the pseudorange and the rate of change of the pseudorange, 
said matched filter algorithm including: 

(i) demodulating said digitized signal relative to the data sequence. 
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32. The method of claim 3 1 , further including the step of: 
(d) aligning the bits of the data sequence with said digitized signal, prior to 
said demodulation. 


33. The method of claim 32, wherein said aligning is within about one 
millisecond. 


34. The method of claim 3 1 , wherein said matched filter algorithm further 
includes: 

(ii) arranging said digitized signal as columns of a matrix that includes a 
plurality of rows, each said column including consecutive said bits of 
said digitized signal that correspond to an integral number of frames of 
the pseudonoise sequence; and 

(iii) performing a discrete orthogonal transform on each said row. 


35. The method of claim 34, wherein said matrix includes a plurality of 
elements, and wherein said matched filter algorithm further includes: 

(iv) multiplying each said element by a respective Doppler compensation 
factor. 


36. 


includes: 


(v) 


The method of claim 35, wherein said matched filter algorithm further 
convolving each said column with the pseudonoise sequence; and 
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(vi) identifying at least one peak in said matrix, each said at least one peak 
having a row coordinate and a column coordinate. 

37. The method of claim 36, wherein said receiving, said digitizing, said 
arranging, said performing, said multiplying, said convolving and said identifying are 
effected on a plurality of instances of the transmitted signal, and wherein said matched 
filter algorithm further includes: 

(vii) inferring the pseudorange and the rate of change of the pseudorange 
from said row coordinates and said column coordinates. 

38. The method of claim 36, wherein said receiving, said digitizing, said 
arranging, said performing, said multiplying and said convolving are effected on a 
plurality of instances of the transmitted signal to produce a corresponding plurality of 
said matrices, and wherein said matched filter algorithm further includes: 

(viii) integrating said matrices non-coherently to produce a sum matrix; 

(ix) identifying at least one peak in said sum matrix, each said at least one 
peak having a row coordinate and a column coordinate; and 

(x) selecting, from among said at least one peak, a most likely peak, said 
row coordinate of said most likely peak then corresponding to the 
pseudorange, and said column coordinate of said most likely peak then 
corresponding to the rate of change of the pseudorange. 
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ABSTRACT OF THE DISCLOSURE 
A method and system for finding the position of a mobile unit with respect to 
the satellites of a satellite network such as the Global Positioning System and with 
respect to the base stations of a wireless communications network. Each satellite 
transmits a signal that consists of a series of frames of a pseudonoise sequence. The 
frames of a signal received from the satellite network by the mobile unit are arranged 
as columns of a matrix and are processed coherently to provide estimated 
pseudoranges and estimated rates of change of pseudoranges for in view satellites. 
The coherent processing includes performing an orthogonal transform on the rows of 
the matrix, multiplying the elements of the matrix by Doppler compensation factors, 
and then, for each satellite in view, convolving the columns of the matrix with the 
pseudonoise sequence of that satellite. Other pseudoranges are inferred from 
synchronization burst sequences received by the mobile unit from one or more base 
stations. If the base stations are not mutually synchronized, reference units are 
provided to determine the time offsets between the base station clocks and the satellite 
network clock and to transmit these time offsets to the mobile unit, so that all the 
pseudoranges include the same bias relative to the satellite network clock. The 
pseudoranges are processed and reconciled by a dynamic estimator that includes a 
linear optimal estimator such as a Kalman filter, to produce successive estimates of 
the location of the mobile unit as the mobile unit moves. 
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